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Abstract 

This paper deals with a general class of integrals, the particular cases of 
which are connected to outstanding problems in astronomy and physics. Re- 
action rate probability integrals in the theory of nuclear reaction rates, Kratzel 
integrals in applied analysis, inverse Gaussian distribution, generalized type-1, 
type-2 and gamma families of distributions in statistical distribution theory, 
Tsallis statistics and Beck-Cohen superstatistics in statistical mechanics and 
the general pathway model are all shown to be connected to the integral under 
consideration. Representations of the integral in terms of generalized special 
functions such as Meijer's G-function and Fox's H-function are also pointed 
out. 

Keywords: H-function, Kratzel integral, generalized beta, gamma and inverse 
Gaussian densities, pathway model. 

1. Introduction 

In this paper we will consider a general class of integrals connected with 
the pathway model of Mathai (2005). These will enable us to address a wide 
class of problems in different areas such as inverse Gaussian processes in the 
area of stochastic processes, Kratzel integral in applied analysis, generalized 
type-1, type-2 and gamma densities in statistical distribution theory, Tsallis 
statistics in non-extensive statistical mechanics, superstatistics in astrophysics 
problems, reaction probability integrals in nuclear reaction rate theory, and 
many related problems, which may be seen from the formalism introduced in 
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this paper. Consider the following integral: 

poo 

f(z 2 \ Zl )= / x" / - 1 [l + z 5 1 (a-l)x s }-^[l + z p 2 (/3-l)x- p }-^ 1 (1) 
Jo 

for a > 1, & > 1, z x > 0, z 2 > 0, 5 > 0, p > 0, £(7 + 1) > 0, 

1 T + 1 1 1 

R ^— r " > °' R ttt " > 

a — 1 p — I p 

poo 1 ~ 

= / -/i(x)/ 2 (^)dx (2) 

J X X 

where 3ft (•) denotes the real part of (•). 

h{x) = x*[l + z[{a - l)/]-^r, / 2 (x) = [1 + (/? - lKn^ 1 (3) 
with Mellin transforms 

M h { S ) = [6zr s (a-i)^r 5 j f !l_\ - * w 

SR( 7 + s) > 0, 3? (— L. - li-f") > 
\a — 1 / 

and 



n^rc— — -) 

1 V/3-lJ 

&(s) > 0, 3ft f -J— >0. 

Hence the Mellin transform of f(z 2 \ Zl ), as a function of z 2 , with parameter s 
is the following: 

M f ( Z2 \ Zl) (s) = M h {s)M h {s) 

dz^a-l)^ r(^) 

for &( 7 + s) > 0, '— - > 0, &(s) > 0, 

a — 1 

3ft(-^---) >0,*i >0,z 2 >0. 
P - 1 P 
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Putting y — i in (1) we have 

r OO — -y 

Mh) = / — [1 + z((a - + - l)j/T^dy. (7) 

Jo y 

Evaluating the Mellin transform of (7) with parameter s and treating it as a 
function of zi, we have exactly the same expression in (6). Hence 

MfWzM = M fMza )(s) = right side in (6). (8) 

By taking the inverse Mellin transform of Mf( Z2 \ Zl )(s) one can get the integral 
f(z 2 \zi) as an H-function as follows: 
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f(z 2 \z 1 ) = c 1 H 22 
where 



Zl z 2 (a-1) H(3-l 



c = 5p Z J(a — 1) « . 



K 1 a-l+l'J)^ 1 /9-1'p) 



(9) 



where H™^ n (-) is a H-function which is defined as the following Mellin-Barnes 
integral: 

I {a\,a-i),...,(a p ,a p ) 
l(bi,/3i),...,(6 9 A) 



Tjm,n 



where 



= ^/^(-)d- do) 

, n _ {n j "L 1 r(& J + /3^)}{n; =1 r(i-q,-M} nn 

{n J i m+ i r ( 1 - 6 J -ft s )}{nUi%+^ s )} 

where L is a suitable contour, = l,...p,/3j,j = 1, ...,<? are real positive 
numbers, = 1, ...,q,a,j,j = l,...,p are complex numbers and L separates 
the poles of F(bj + f3js),j = 1, ...m from those of T(l — aj — a.js),j = 1, n. 
For more details about the theory and applications of H-function see Mathai 
and Saxena (1978) and Mathai et al. (2009). 

The integral in (1) is connected to reaction rate probability integral in 
nuclear reaction rate theory in the non-resonant case, Tsallis statistics in 
non-extensive statistical mechanics, superstatistics in astrophysics, general- 
ized type-2, type-1 beta and gamma families of densities and the density of a 
product of two real positive random variables in statistical literature, Kratzel 
integrals in applied analysis, inverse Gaussian distribution in stochastic pro- 
cesses and the like. Special cases include a wide range of functions appearing 
in different disciplines. 

Observe that fi(x) and f 2 (x) in (3), multiplied by the appropriate normal- 
izing constants can produce statistical densities. Further, fi(x) and f 2 (x) are 
defined for — oo < a < oo, — oo < f3 < oo. When a > 1 and z\ > 0, S > 0, fi(x) 
multiplied by the normalizing constant stays in the generalized type-2 beta 
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family. When a < 1, writing a — 1 = —(1 — a), a < 1 the function fi(x) 
switches into a generalized type-1 beta family and when a — > 1, 



lim ^(x) = e"*** (12) 



a— ¥l 



and hence fi(x) goes into a generalized gamma family. Similar is the behavior 
of f 2 (x) when (3 ranges from — oo to oo. Thus the parameters a and j3 cre- 
ate pathways to switch into different functional forms or different families of 
functions. Hence we will call a and j3 pathway parameters in this case. Let us 
look into some interesting special cases. Take the special case /3 — > 1, 

POO 

f 1 (z 2 \z 1 )= / a; 7 - 1 [l + ^(a-l)a:*]-=^e-^ :e "''da; (13) 
Jo 

a > 1, z 1 > 0, z 2 > 0, S > 0, p > 0. Put y = \ 

poo 

fM\z 2 ) = / y-^l + 4(a - l)y- s ]-^e-^ P dy (14) 
a > 1, zi > 0, z 2 > 0, 8 > 0, p > 0. Let a ->■ 1 in (1) 

/•oo 

/ 2 (^ki)= / a^-V^l + ^Ctf- l)ar-"]"^da; (15) 

JO 

/5 > 1, zi > 0, z 2 > 0, 5 > 0, p > 0. 

roo 

f 2 ( Zl \z 2 ) = / x-T-^-^'fl + 4(0 - l)x p }-—idx (16) 
Jo 

P > 1, > 0, z 2 > 0, 5 > 0, p > 0. Take a ->■ 1, /? ->■ 1 in (1) 

POO 

f 3 (z 2 \ Zl ) = / x^e-^-^-'dx (17) 
•/ o 

> 0, z 2 > 0, 5 > 0, p > 0. 

POO 

f 3 ( Zl \z 2 )= / x^-V^'-^dx (18) 
Jo 

zi > 0,z 2 > 0,5 > 0,p > 0. 



In all the integrals considered so far, we had one pathway factor containing 
x s and another pathway factor containing x~ p , where both the parameters 
5 > and p > 0, in the integrand. Also the integrand consisted of non- 
negative integrable functions and hence one could make statistical densities out 
of them. In statistical terms, all the integrals discussed so far will correspond 
to the density of u — X\X 2 where X\ and x 2 are real scalar random variables, 
which are statistically independently distributed. 

Now we will consider a class of integrals where the integrand consists of 
two pathway factors where both contain powers of x of the form x s and x p 
with both 5 and p positive. Such integrals will lead to integrals of the following 
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forms in the limits when the pathway parameters a and ft go to 1: 



x 7 e 



a > 0, b > 0, 5 > 0, p > 0. Observe that the evaluation of such an inte- 
gral provides a method of evaluating Laplace transform of generalized gamma 
densities by taking one of the exponents 5 or p as unity. Consider the integral 

POD 

h= x 7 [l + z((a-l)x s ]-^[l + z p 2 {ft -l)x p r~^dx, (19) 
Jo 

a > I, ft > l,zi > 0,z 2 > 0,6 > 0,p > 0. Since the integrand consists of 
positive integrable functions, from a statistical point of view, the integral I± 
can be looked upon as the density of u — ^, where X\ and X2 are real scalar 
random variables which are independently distributed or it can be looked upon 
as a convolution integral of the type 

poo 

/ vh(uv)f 2 {v)dv (20) 
Jo 

Let us take 

„si-— 



fi(xi) = ci[l + (a - l)^] °-i,u = z 1 
/ 2 ( X2 ) = ^^[l + z p 2 {ft - l)a%-T& 



Taking the Mellin transforms and writing as expected values 



E( Xl ) s - 1 = ci / x{~ *[1 + (a - l)x d 1 ]"^dxi 



./0 




Cl 




<5(a — 


l)f 


POD 




02 / 


7— 


./0 





-,ft(s) > 0,SR( 



a — 1 5 



>o 



c 2 



7— s+1 j p / 1 7— s+1 

p ) U-i p 



■Y — S + 1 

p[^(/3-l)] — 



(A 



5R(7 - s + 1) > 0, 5ft(-^r " 7 S + 1 ) > 0. 

ft -I p 
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Therefore the density of u = 21 is given by 

X2 



g(u) = 



ClC 2 



p ( 7+1 _ A r ( _j 7+1 _i_ A 

1 f \ p p) p p) 



Therefore 



5p[z p 2 ((3-l)]- 
1 



27T2 . /l 



ClC 2 



r (A) r (^) 

z 2 (ck - 1)« 



M^-i)]^r(^)r(^) 



x # 2 ' 2 

x ^2,2 



z 2 (a - 1)> 



MP- 1 )* 



i 1 1+ P 'D'i 1 a-l'l) 

Co M ^ ^ 



Z" 00 i 
/ 4 = / x 7 [l + 4(a-l)^]"^[l + 4(/3-l)^]"~da;, 
Jo 



a > l,/3 > 1,5 > 0,p > 
1 

M^(/3-i)]^r(^)r(^ 



X 11 2,2 



z 2 {a - l)s 



(°'i)'(^r~ 2 7 i 'p) 



MP- 1 )' 

Now by putting y = | we can get an associated integral 



✓•oo 

J 4 = / y-^tl + ^a-l^-'l-^rfl + ^-l^-^dj/. 

./0 



Now, we can look at various special cases of lim^! or lim^! or lim^i^i. 
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These lead to some interesting special cases. 



hi — lim J 4 



'4.2 



/ 

JO 



x J e~ z * x [1 + z p 2 ((5 - l)x p r— dx 



lim I4 



rr2,l 
^1,2 



^2 



7+1 i % 



(°'J)'(^T _2 7 i ')) 



/ 

Jo 



x 7 [l + 4 (a - l)x 6 Y^e~ z * xP &x 



1.2 



h.3 = 



p^ +1 r (^) 

lim ^4 

a-*-l,/3->-l 



#2,1 



2; 2 (a - l)<s 



(04) 



7g-(^ia:)' 5 -(^22:) P c J x 



P<^2 



7+1 



^2 









- 2 e 4 x s 4 x P ci x . 



When a < 1 and j3 < 1 also we can obtain corresponding integrals, which are 
finite range integrals, by going through parallel procedure. In this case the 
limit of integration will be < x < max{ei, e 2 } where e\ = [zf(a — and 
e 2 = \z p 2 {fi — 1)]"- The details of the integrals will not be listed here in order 
to save space. 

1.2. Case of a < 1, or f3 < 1. 

When a < 1, writing a — 1 = — (1 — a) we can define the function 

9l ( x ) = x i[i + z [{a - l)x s ]-^ =x 7 [l - zf(l- a)x 5 ]^,a < 1 (21) 



for [1 — z((l — a)x 5 ] > 0,a < 1 =>- x < — r and g±(x) = elsewhere. In 

zi(l— a)S 
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this case the Mellin transform of gi(x) is the following: 



p<X> r- y 

h x (s) = / x s - 1 g 1 (x)dx = / Zl(1 - Q)7 xi+'-^l - 4(1 - a)x 5 ]^dx (22) 
Jo Jo 

i r(^)r(-^- + i) 

= 1 , , ' ,fth + s)>0,a<l,5>0. (23) 

Then the Mellin transform of f(z 2 \zx) for a < l,/3 > 1 is given by 

rfe + i) r(2+2) r(j)r(^-j) 



M* a |*i( 5 ) 



<^S*7 + '(/? - i)5(i - «)^ r(¥ + ^ + 1) r(^ 



13-1) 



K( 7 + s) > 0,ft(s) > 0,3? 



/5- 1 P, 



Hence the inverse Mellin transform for a < 1 , f3 > 1 is 



^7(i-a)ir(^) 



X #. 



2,1 
2,2 



lim /(zsl^i) = ^ ° + 7x H l? 



lim f(z 2 \zi) 



1 



a— >1 



P^(^)^ 



zi2 2 (l - a) s (/3 - l) p 

2:1^2(1 - a 

W/? " 1) 



(o^)>(la) 



2,1 



#1,2 



i|(l+T^ + ?.j) 



(o,i),a,i) 



lim f(z 2 \ Zl ) 



a-H,P->l 



pSz{ 



2,0 
7 "0,2 



ZlZ 2 \ 



(o,±),(|,f) 



(24) 



(25) 
(26) 
(27) 
(28) 



In /(^l^i) if P < 1 we may write (3 — 1 = —(1 — /3), and if we assume 
[1 — ^(1 — /3)x~ p ] 1 ^ > =>- x > z 2 (l — f3)p then also the corresponding 
integrals can be evaluated as H-functions. But if a < 1 and j3 < 1 then from 
the conditions 

l-z{{l-a)x & > =>• x < r and \-z p 2 {\-p)x~ p > =>• x > z 2 (l-p)p 

zi(l — a)<s 

and the resulting integral may be zero. Hence, except this case of a < 1 and 
j3 < 1 all other cases: a > 1, /? > 1; a < 1, /3 > 1; a > 1, /? < 1 can be given 
meaningful interpretations as H-functions. Further, all these situations can 
be connected to practical problems. A few such practical situations will be 
considered next. 
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2. Specific Applications 

2.1. Kratzel Integral 

For 5 = 1, z£ = z, Z\ = 1 in fz(z 2 \zi) gives the Kratzel integral 

poo 

f 3 (z 2 \ Zl ) = / x^e^-^dx (29) 
Jo 

which was studied in detail by Kratzel (1979). Hence / 3 can be considered 
as generalization of Kratzel integral. An additional property that can be seen 
from Kratzel integral as fa is that it can be written as a H-function of the 
type Hq^(-). Hence all the properties of H-function can now be made use of 
to study this integral further. 

2.2. Inverse Gaussian Density in Statistics 

Inverse Gaussian density is a popular density, which is used in many disci- 
plines including stochastic processes. One form of the density is the following 
(Mathai, 1993, page 33): 

f( x ) = c x'h~^ + ^,fi ^ 0, x > 0, A > (30) 

i >> 

where c = ir~2e^ . Comparing this with our case f3{Zi\z 2 ) we see that the 
inverse Gaussian density is the integrand in fz(zi\z 2 ) for 7 = |, p = 1, z 2 = 
= I'^i = f- Hence / 3 can be used directly to evaluate the moments 
or Mellin transform in inverse Gaussian density. 

2.3. Reaction Rate Probability Integral in Astrophysics 

In a series of papers Haubold and Mathai studied modifications to Maxwell- 
Boltzmann theory of reaction rates, a summary is given in Mathai and Haubold 
(1988) and Mathai and Haubold (2008). The basic reaction rate probability 
integral that appears there is the following: 

h= x^ x ^ ax - zx ^ dx. (31) 
Jo 

This is the case in the non-resonant case of nuclear reactions. Compare integral 

1\ with jj,{z 2 \z^). The reaction rate probability integral l\ is fz{z 2 \zi) for 
1 

S = 1, p = |, z 2 = z. The basic integral ii is generalized in many different 
forms for various situations of resonant and non-resonant cases of reactions, 
depletion of high energy tail, cut off of the high energy tail and so on. Dozens 
of published papers are there in this area. 

2.4. Tsallis Statistics and Superstatistics 
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Tsallis statistics is of the following form: 

/ x (a;) = c 1 [l + (a-l)a;]-^T. (32) 

Compare f x (x) with the integrand in (1). For z 2 = 0, 5=1,7 = 1 the 
integrand in (1) agrees with Tsallis statistics f x {x) given above. The three 
different forms of Tsallis statistics are available from f x (x) for a > l,a < 
l,a — > 1. The starting paper in non-extensive statistical mechanics may be 
seen from Tsallis (2009). But the integrand in (1) with z 2 — 0, Z\ — 1, a > 1 is 
the superstatistics of Beck and Cohen, see for example Beck and Cohen (2003), 
Beck (2006). In statistical language, this superstatistics is the posterior density 
in a generalized gamma case when the scale parameter has a prior density 
belonging to the same class of generalized gamma density. 

2.5. Pathway Model 

Mathai (2005) considered a rectangular matrix-variate function in the real 
case from where one can obtain almost all matrix-variate densities in current 
use in statistical and other disciplines. The corresponding version when the 
elements are in the complex domain is given in Mathai and Provost (2006). 
For the real scalar case the function is of the following form: 

f(x) = c*\x\' 1 [l-a{l-a)\x\ 5 ]^ (33) 

for — oo <x<oo,a>0,r]>0,5>0 and c* is the normalizing constant. 
Here f(x) for a < 1 stays in the generalized type-1 beta family when [1 — 
a(l — a)|x| 5 ]i^ > 0. When a > 1 the function switches into a generalized 
type-2 beta family and when a — > 1 it goes into a generalized gamma family 
of functions. Here a behaves as a pathway parameter and hence the model is 
called a pathway model. Observe that the integrand in (1) is a product of two 
such pathway functions so that the corresponding integral is more versatile 
than a pathway model. Thus for z 2 = in (1) the integrand produces the 
pathway model of Mathai (2005). 
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